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Abstract. Some recent observations of the abundances of s-process, r-process, and a elements in metal-poor stars 
have led to a new scenario for their formation. According to this scenario, these stars were born in a globular cluster 
and accreted the s-process enriched gas expelled by cluster stars of higher-mass, thereby modifying their surface 
abundances. Later on, these polluted stars evaporated from the globular cluster to constitute an important fraction 
of the current halo population. In addition, there are now many direct observations of abundance anomalies not 
only in globular cluster giant stars but also in subgiant and main-sequence stars. Accretion provides again a 
plausible explanation for (at least some of) these peculiarities. Here we investigate further the efficiency of the 
accretion scenario. We find that in concentrated clusters with large escape velocities, accretion is very efficient 
and can indeed lead to major modifications of the stellar surface abundances. 
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1. Introduction 

Strong correlations between the r-process and s-process el- 
ement abundances and the a-element abundances in field 



metal-poor stars have been reported (Jchin ct al. 1998, 
1999), separating these stars into two sub-populations. 
The r-process elements correlate linearly with the a- 
elements, with a clumping at the maximum value of 
[a/Fe\. The s-process elements, on the other hand, ex- 
hibit a more complex behavior when plotted against the 
a elements, and form a "two-branches diagram" , which 
we show schematically in Fig. 1. We emphasize here that 
the observed stars being dwarf stars, they cannot have 
synthetized these s-elements in their interior. 

In order to explain these observ ations the EASE 
scenario has been developed ( Ichin ct al. 1999 ; 



Parmcnticr ct al. 1999| , ^0004 |2000bD which lii^s 



PopIIb 



EASE 



Popllfc 



the metal-poor field halo stars to the halo globular 
clusters (GCs). According to this model, the field halo 
stars are born in GCs. The evolution of the GCs is 
separated into two phases, directly connected to the 
two subpopulations now observed in the halo. During 
the first phase, "first generation" stars are formed. The 
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Fig. 1. The "two-branches diagram" (see |Jehin et al. 1999 
for details). 



initial mass function in the very metal-poor medium 
favors the high masses, with 3 Mq < M < 60 M©. The 
most massive of these stars quickly evolve to become 
supernovae, ejecting a- and r-process elements into the 
surrounding intracluster medium. The explosions trigger 
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the formation of an outwardly expanding supershell. 
The primordial ISM gas is swept into the shell with the 
supernovae ejecta, and the resulting gas is enriched in a- 
and r-process elements. A new burst of star formation 
occurs in the shell, where the density is very high. This 
triggered star formation presumably has a Salpeter-like 
initial mass function. A proto-globular cluster is born. If 
the gravitational potential is too weak, it gets disrupted 
before having the chance of forming a GC, and the stars 
become field halo stars of the PopIIa, as defined on Fig. |l|. 
The metallicity reached in the proto-GC depends on the 
number of supernovae which have exploded before the 
proto-GC gets disrupted. The a-element abundances are 
fixed by that time as well as by the initial mass function 
(supernova progenitors of different masses produce differ- 
ent relative amounts of a elements). The proto-globular 
cluster may in some cases survive this supernova phase. 
If the stars fall back under the effect of gravity, they 
form a GC. Detailed dynamical ca lculations have been 
performed ( Parmcntier et al. 199£ ) which show that 
proto-globular clusters can indeed survive the supernova 
phase and metallicities such as those observed today in 
GCs can be achieved in this scenario. It has also been 
shown that the EASE scenario, with its two separate 
star formation events, does not contradict observational 
constraints such as the narrow metallicity range found in 
any given cluster. 

During the second phase, intermediate-mass (sec- 
ond generation) stars evolve, becoming asymptotic gi- 
ant branch (AGB) stars. Those stars produce s-elements, 
which are brought to the surface through the so-called 
third dredgc-up occurring after the thermal pulses. The s- 
elements are then ejected into the surrounding intracluster 
medium by stellar winds. The intracluster gas is therefore 
enriched in s-elements, while the a- and r-process element 
abundances have been fixed at their highest value after all 
massive stars have exploded. If the cluster can retain this 
enriched gas, it can be accreted by the lower-mass stars 
which are still on the main sequence. If the accreted ma- 
terial does not get mixed within the entire star but only 
in its convective zone, which contains about 1% of the 
total stellar mass for stars still near the main-sequence 
turn-off of present GCs, then even a small amount of ac- 
cretion can lead to an appreciable enrichment of the stel- 
lar atmosphere. When the star evolves towards the giant 
stage, its convective zone becomes much larger, and the 
s-element enrichment is diluted considerably. We point 
out here that some mechanisms could be at work to al- 
ready trigger the dilution of the accreted matter below 
the convective zone of main- sequence stars (Proffitt 1989; 
Proffitt & Michaud 1989), thus diminishing the impact on 
the surface abundances of the deep convective envelope de- 
veloping in red giant stars. We come back to this issue in 
more detail in Sect. 2. The accreting star can be ejected 
out of the GC at any time, by evaporation or by tidal 
forces possibly leading to the total disruption of the clus- 
ter while crossing the disk. There is indeed now observa- 
tional evidence that GCs lose stars: two cannon-ball stars 



in 47 Tuc ([Meylan fc al. 199l[ [Jorissen fc Meylan 2001 



see also Capaccioli et al. 1993; Piotto et al. 1997, and ref- 
erences therein), and, above all, the recent observations 
of "tidal tails" around many GCs (Grillmair et al. 1995; 



Leon ct al. 2000| ; pdcnkirchcn et al. 2001^ Stars 



hanced in s-proccss elements by accretion and then 
stripped from GCs form the Pop lib halo. 

The hypothesis that the gas ejected by the 
intermediate-mass stars in GCs can be accreted by 
the other (lower-mass) cluster stars has been ex- 
amined before. There is strong observational and 
theoretical evidence that stars reaching the AGB 
phase lose a large amount of mass. Several studies 
have been devoted to the fate of this gas in clus- 
ters 



VandenBert 



(IScott fc Rose 1975 



Faulkner 1977 



Scott & Durisen 1978 



Faulkner & Freeman 1977 




VandenBerg 1978 






Faulkner 1984 



Faulkner fc Coleman 1984 ; Smith 1996 ). Stellar ejecta in 
GCs with shallow potential wells can leave the cluster via 
a smooth wind-like outflow. Faulkner fc Freeman (1977) 
and VandenBerg fc Faulkner (1977) find steady-state 
time-independent fiow solutions in clusters of 10^ M©. 
However, they show that in GCs with deep enough 
potential wells, the gas can accumulate into the cluster, 
forming a central reservoir with a radius comparable to 
the GC core radius. Faulkner fc Coleman (1984) show 
that a small number of low-velocity low-mass stars in 
the cluster core can in this case accrete enough matter 
to form 10 M0 black holes, which can be considered 
as an extreme case. All those studies looked at what 
happens in present-day GCs. Smith (1996) addressed the 
question whether such winds might have been possible 
within young GCs, during epochs of much higher stellar 
mass-loss rates. Many GCs have high enough escape 
velocities to retain at least some of the stellar ejecta. 
A substantial amount of intracluster material could 
have been acquired when the turn-off mass was about 
5M0. In less tightly-bound clusters, the stellar ejecta 
are lost from the cluster either stochastically or through 
a continuous wind. In the cases where the stellar ejecta 
are retained in the cluster, the gas could be accreted 
by other cluster stars, thereby modifying their surface 
composition. Several models estimating the efficiency of 
the accretion scenario have appeared in the literature (see 
e.g. D'Antona ct al. 198S; Smith 1996). A quantitative 



model investigating the recycling of nova ejecta has 
been presented by Smith & Kraft (1996). Qualitative 
discussions of the accretion model have been presented by 
Bell et al. (1981), Norris fc Da Costa (1995), and Cannon 
et al. (1998). 

In this paper, we study the process of accretion onto 
GC stars from a central reservoir of gas. In particular we 
take into account the evolution of this process with time. 
Indeed the rate of mass ejection into the cluster is strongly 
time-dependent, as it is highest when intermediate-mass 
stars reach the mass-losing stage, decreasing very rapidly 
as the turn-off mass decreases. In Sect. 2, we examine the 
possible links between halo stars and GCs. We address the 
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problem of the "missing" intracluster gas in Sect. 3. We 
present some observational signatures that may possibly 
be related to accretion by cluster stars in Sect. 4. In Sect. 5 
we present the method used to calculate the efficiency of 
accretion of gas by GC stars. The results are summarized 
in Sect. 6. Finally Sect. 7 contains the conclusions. 

2. Unease with the EASE scenario 

Although the EASE scenario nicely explains some abun- 
dance correlations in halo stars through a sequence of 
events in GCs, it faces several problems. 

The GC evaporation is not the only way to populate 
the halo. In particular, the discovery that the Sagittarius 
dwarf galaxy is in the process of being tidally disrupted 
by our Galaxy and its stars incorporated in the halo has 
led to a serious reconsideration of the role of mergers in 
the formation of the halo. Initially put forward by Searle 
& Zinn (1978), in opposition to the classical rapid col- 
lapse envisaged by Eggen et al. (1962), the merger sce- 
nario has gradually gained consensus. Relics of these past 
mergers may be the streamers highlighted by Majewski 
et al. (1996) or the blue metal-poor stars highlighted by 
Preston et al. (1994). 

If field halo stars come from globulars, they should ex- 
hibit the same pattern of abundance peculiarities as GG 
stars [like C-(Na,Al) anticorrelation, N-Na-Al correlation]. 
Hanson et al. (1998) have shown that, based on their lo- 
cation in the ([Na/Fe], [Mg/Fe]) diagram, halo field and 
GG stars of the same metallicity are not surrogates of one 
another. In particular, no field halo stars seem to exhibit 
the large [Na/Fe] overabundances observed in some bright 
giants in GGs (Pilachowski et al. 1996). Nevertheless, the 
scatter in [Na/Fe] a t a given metallicity observed among 
halo stars (Fig. 4 of Hanson ct al. 1998 ) is certainly large 
enough to accomodate the 0.3 to 0.4 dex range observed 
in 47 Tuc main-sequence stars and attributed to the wind 
pollution scenario. 

If the accretion scenario is at work in GGs, one would 
expect that the pollution levels decrease from main se- 
quence to giant stars as the convective envelope deepens, 
and dilutes the material initially restricted to a thin sur- 
face layer. Several among the chemical anomalies observed 
in GGs and mentioned above do not exhibit this trend, 
however, as they remain unchanged all the way from the 
main sequence to the giant branch. This suggests that the 
accreted matter is rapidly diluted even in main sequence 
stars. At least two processes may be invoked to trigger 
the dilution of the accreted matter into the deep radia- 
tive layers of main sequence stars. A turbulent diffusion 
mechanism triggered by the inversion of mean molecular 
weight due to the higher He content of the accreted mat- 
ter has been presented by Proffitt (1989) and Proffitt & 
Michaud (1989) in relation with Am and Ba stars. These 
authors show that even a modest composition inversion 
(say AX = -0.001, and correspondingly AY = +0.001, 
where X and Y are the hydrogen and helium mass frac- 
tions, respectively), over the outer 10"'^ Mq of the star 



results in dilution of the added layer by a factor of 50 in 
10® y. The depth reached by such a mixing may be suffi- 
cient in order to prevent further dilution during the first 
dredge-up. This effect may explain why the level of abun- 
dance anomalies is similar in dwarf Ba stars and in giant 
Ba stars (North et al. 1994), and why the level of abun- 
dance anomalies does not change between dwarf and giant 
stars in GGs. Another cause of mixing specific to GGs is 
tidal mixing caused by occasional close encounters of stars 
in clusters as suggested by Gox (1998). 

The present paper is motivated by the connection be- 
tween s-process enhancements in field halo stars and in 
cluster stars that evaporated from the cluster. There is 
not much evidence, however, for s-process anomalies due 
to accretion in globular cluster stars, except for uj Gen 
(but see the discussions of Sects. 4.1 and ^ and for the 
metal-poor cluster M15 where Sneden et al. (1997) re- 
ported bimodal abundance distributions for both [Ba/Fe] 
and [Eu/Fe] (with a strong correlation between Ba and 
Eu) in the bright giants. 



3. Low level of intracluster gas 

The absence of interstellar material in GGs ha s 
been a long-st anding puzzle (e.g. Roberts 1988 ). 



Freire et al. 2001 recently achieved a positive detec- 
tion of 0.1 M0 of ionized gas in the core of 47 Tuc. 
But this mass is much less than the ~ 100 Mq of 
intracluster material expected to accumulate between two 
s uccessive pass ages of the GG through the galactic disk 
( [Roberts 1988| ). Knapp et al. (1995, 1996) have recently 
reinvestigated this question by setting upper limits on 
both the dust and ionized gas contents for a number of 
GGs using IRAS and radio observations, respectively. 
Knapp et al. (1996) conclude that their 'upper limits on 
the mass of ionized gas are, for all observed clusters, much 
smaller than the amount expected to have accumulated 
from stellar mass loss if all that gas remains in the cluster 
between galactic plane passages'. Simple models of the 
equilibrium velocity distribution of the gas supposedly 
photoionized by post-AGB stars indicate, however, that 
this gas is able to escape from the cluster in the time sep- 
arating two passages through the galactic plane. Taking 
this effect into account, the amount of gas predicted to 
accumulate from stellar mass loss in this time is now well 
below the observational upper limits. As far as the dust is 
concerned, Knapp et al. (1995) conclude that the upper 
limits are much lower than the amounts expected to have 
accumulated over 10® y from stellar mass loss. 

Recycling of the mass lost by the evolved stars, by 
the accretion scenario investigated in the present paper, 
offers another way (see e.g., Freire ct al. 2001 for a list of 



some of the other possible mechanisms) to account for the 
apparent lack of intracluster matter. 



4 



A. Thoul et al.: Accretion of gas by globular cluster stars 



4. Abundance anomalies observed in globular 
cluster stars, and their possible relevance to 
the accretion scenario 

4.1. Peculiar red stars in lu Cen 



The GC uj Cen is discussed separately as it bears sev- 
eral peculiarities. It is the most massive and largest GC 
( Harris 199(\ ) and exhibits a large spread in metallic- 
ity (-1.8 < [Fe/H] < -0.8; [Norris fc Da Costa 1995"| ) 
with some spatial asymmetry (Jurcsik 1998), but the 
most significant peculiarity of lo Cen in the present 
context i s undoubtedly the p resence of CH, B a and 
S sta rs ([Dickens fc Bell 19761 [Lloyd Evans 1974 1983, 



1986; [Vanture et al. 1994[) , as listed in Table [l[. Mild 



overabundances of the heavy elements produced by the 
s-process of nucleosynthesis are observed in the pe- 
cuHar red stars (PRS) of u Cen (Table |l|). With 
the exception of the CH star ROA 279, PRS in 
oj Cen form the tail of the general trend of increas- 



typical of lu Cen ( 


Vanture et al. 1994 


; Smith ct al. 1995; 


Norris fc Da Costa 1995 


). PRS in 10 Cen should thus not 



be considered as exceptional, but rather as the most ex- 
treme members of a chemically-inhomogeneous popula- 
tion. The PRS in cu Cen generally turn out to be some- 
what enriched in carbon, since the u Cen giants studied 

by Norris & DaCosta (1995) typically have [C/Fe] 0.6, 

as compared to [C/Fe] values as large as -1-0.2 for some of 
the PRS listed in Table |^. The existence of such carbon- 
enriched giants is another specific feature of uj Cen, since 
in GGs with [Fe/H] ~ —1.0, red giants within 1 mag of 
the red giant tip exhibit carbon depletions by factors in 



the range 0.3-1.0 dex instead (Briley et al. 1992, and ref- 
erences therein). 

The PRS of u! Cen are fainter than the luminosity 
threshold of the thermally-pulsing AGB {M^ ^ —2.5), 
according to the M = 1.25 M0,Z = 0.001 evolutionary 
track of Charbonnel et al. (1996) [their Fig. 2]. Therefore, 
internal nucleosynthesis on the TP- AGB cannot account 
for their chemical peculiarities. 

Barium and CH stars in the field have been 
shown to belong systematically to binary systems 
( McClure fc Woodsworth 1990 ), with their chemical 
anomalies resulting from the accretion of carbon- and 
s-process-rich matter from their AGB companion that 
has now evolved into a white dwarf. A radial-velocity 
monitoring of the PRS in lu Cen has been performed 
by Mayor et al. (1998), and has revealed that with the 
exception of the classical CH stars ROA 55 and ROA 70, 
as well as possibly the S stars ROA 320 and V6, PRS 
in OJ Cen do NOT belong to binary systems, however 
(Mayor et al. 1996; Table A scenario specific to GGs 
seems thus at work to produce these PRS. One possible 
scenario is that those PRS were formerly member of 
binary systems that were disrupted by a close encounter 
with another cluster star. Mayor et al. 199(^ estimate that 
systems with periods in excess of 3 to 10 y, depending 



on their position in the cluster, have been disrupted by 
tidal encounters over the lifetime of the cluster. Before 
being disrupted, Ba, S or CH systems with such long 
periods may indeed have been produced by the same 
mass-transfer scenario as that invoked to account for the 



PRS observed in the field ( [Jorissen et al. 1998[ 1999). 

Norris fc DaCosta (1995) and Gannon et al. (1998) 
have suggested instead that the PRS envelope has been 
polluted by winds expelled in the intracluster medium by 
an early generation of AGB stars. However, the simulta- 
neous increase of Fe and the s-process elements observed 
in LU Cen requires as well a source of Fe that cannot be 
delivered by AGB stars, a statement not easy to reconcile 
with the pollution scenario, as discussed by Smith et al. 
(1995). Such peculiarities in lu Gen are more in favor of 
a merging between two different clusters, as further dis- 
cussed in Sect. |[ which also shows that accretion is not 
at all efficient in lu Cen. 



4.2. Abundance anomalies in main sequence or 
subgiant stars 

Abundance anomalies have long been known in GC giant 
stars (see e.g. Kraft (1994) and DaCosta (1997) for recent 
reviews). However, the dredge- up processes occurring in 
these stars make it very difficult to disentangle intrinsic 
from extrinsic causes to the observed chemical peculiari- 
ties. Such dredge-up processes are not expected to occur 
in main sequence stars, according to canonical stellar evo- 
lution theories. Abundance anomalies in those stars there- 
fore point in principle to external causes. Because of the 
faintness of main sequence stars in GGs, abundance anal- 
yses of these stars have but only recently become feasible. 
Indications for anomalies have since then rapidly accumu- 
lated. They may be grouped into two classes ~ CN and Na 
anomalies - which are discussed in turn. 



4.2.1. CN anomalies in main sequence stars 

The bimodal nature of the distribution of the CN 
band strength of red giant stars in several GGs has 
been recog nize d since the early 19 80s (see references in 
Kraft 1994 and Cannon et al. 199^ ) . This property has re- 
cently been extended to main sequence stars in 47 Tuc 



( [Bell et al. 1983[ [Briley ct al. 1994[ [Cannon et al. 1998[) , 



M5 and NGG 6752 ([Suntzcff fc Smith 199 1[) , and M13 
and M71 ([Briley fc Cohen 200l[) . The ratio of GN-weak 
to GN-rich stars is similar among dwarf and more lumi- 
nous stars, implying little change in the overall distribu- 
tion of CN with evolutionary state ( Briley et al. 1994 ). 
CN-strong stars tend to be more numerous in the central 
region of 47 Tuc ( Norris fc Freeman 1979 ; Freeman 1985 ). 

Recently, Grundhal et al. (1998) have suggested that 
the significant scatter in the Stromgren ci index observed 
along the red giant branch down to the turnoff for all clus- 
ters more metal-poor than [Fe/H] — -1.2 is attributable 
to star-to-star abundance variations in the GNO elements 
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Table 1. Abundances, luminosities, and binary properties of PRS in uj Cen. Column 1 provides the ROA star 
number from WooUey (1966; last three digits as in WooUey's Tables I and II, or variable numbers from WooUey's 
Table III). Column 2 gives the spectral type from the reference mentioned in column 3, columns 4 and 5 list 
(adopting a distance modulus of 13.92; Harris 1996) and B — V from Mayor et al. 1998, columns 6, 7 and 8 (taken 
from Mayor et al. 1998) give the uncertainty ey^ of the mean radial velocity, the rms value ay^ for stars with more than 
one measurement and the probability -P(x^) that the observed ay^ is due to observational error only. Column 9, labelled 
Vr, provides a flag characterizing the binary nature of the star, according to the following rules: a star is considered as 
spectroscopic binary ('SB') if P{x^) < 0.01 and ay^ > 2.0 km s^^. The latter condition is imposed because red giants 
close to the RGB tip are known to exhibit a radial- velocity jitter of the order of 1 km s~^ (Mayor et al. 1984; see also 
Fig. 1 of Jorissen et al. 1998). If 0.01 < Pix^) < 0.05, or P(x^) < 0.01 and ay^ > 1.5 km s \ the star is considered as 
variable ('var'). Columns 10, 11 and 12 list the [Fe/H], [C/Fe] and s-process abundances when available, with [s/Fe] 
= ([Y/Fc] + [Zr/Fe])/2. Column 13 provides the reference for the abundance data. 



ROA 


Sp. 


Ref 


Mv 


B - V 




Gy 




Vr 


L / J 


fC/Fel 


fs/Fel 

L / J 


Ref 


55 


CH 


1 


-2.43 


1.74 


0.74 


3.98 


0.000 


SB 










70 


CH 


1 


-2.34 


1.82 


0.90 


4.33 


0.000 


SB 










V53 


Ba 


2 


-1.97 


1.78 


0.50 


1.00 


0.015 


var 










201 


S2. 5/1-2 


3 


-1.82 


1.61 


0.52 


1.04 


0.030 


var 


-0.85 


-0.30 


0.33 


4 


219 


Ba 


2 


-1.76 


1.68 


0.41 


1.01 


0.003 


var 


-1.25 


0.22 


0.40 


4 


231 


Ba 


2 


-1.73 


1.63 


0.40 


1.56 


0.000 


var 


-1.10 


-0.60 


0.44 


4 


243 


Ba 


2 


-1.70 


1.31 


0.31 


0.70 


0.155 












248 


Ba 


2 


-1.69 


1.84 


0.46 


1.44 


0.000 


var 


-0.78 




0.17 


4 


270 


Ba 


2 


-1.61 


1.53 


0.37 


0.99 


0.014 


var 


-1.24 


-0.01 


0.34 


4 


276 


Ba 


2 


-1.61 


1.52 


0.37 


0.82 


0.113 












279 


CH 


1 


-1.60 


1.48 


0.33 


0.88 


0.209 




-1.70 




0.48 


4 


286 


Ba 


2 


-1.59 


1.57 


0.25 


0.45 


0.774 












287 


Ba 


2 


-1.59 


1.48 


0.30 


0.67 


0.335 




-1.43 


0.06 


0.26 


4 




S9 /I 9 


O 




1 "iS 


fl 97 


\J.OO 














316 


Ba 


2 


-1.50 


1.66 


0.58 


1.16 


0.017 


var 










320 


S2/1-2 


3 


-1.49 


1.65 


0.86 


2.58 


0.000 


SB? 










321 


Ba 


2 


-1.48 


1.45 


0.23 


0.27 


0.945 












V6 


S3/1 


3 


-1.47 


1.54 


1.18 


2.35 


0.000 


SB? 










324 


Ba 


2 


-1.46 


1.59 


0.25 


0.39 


0.749 












332 


Ba 


2 


-1.45 


1.53 


0.31 


0.75 


0.115 












357 


Ba 


2 


-1.42 


1.48 


0.48 


1.17 


0.000 


var 


-0.85 




0.20 


4 


371 


K5Ba 


3 


-1.38 


1.74 


0.51 


1.25 


0.000 


var 


-0.79 


-0.46 


0.46 


4 




















-1.0 




0.8 


5 


421 


Ba 


2 


-1.25 


1.50 


0.35 


0.86 


0.035 


var 










425 


S3/2 


3 


-1.24 


1.69 


0.33 


0.31 


0.890 












447 


S2/2 


3 


-1.18 


1.71 


0.28 


0.63 


0.338 












451 


Ba 


2 


-1.18 


1.53 


0.41 


0.92 


0.046 


var 










V17 


Ba 


2 


-1.17 


1.80 


0.65 


1.46 


0.000 


var 










480 


Ba 


2 


-1.15 


1.51 


0.49 


0.84 


0.133 




-0.95 


-0.65 


0.38 


4 


505 


Ba 


2 


-1.09 


1.28 


0.27 


0.61 


0.363 












513 


S3/2 


3 


-1.08 


1.67 


0.53 


0.92 


0.065 












577 


CH 


1 


-0.97 


1.55 


0.53 


1.30 


0.008 


var 











References to spectral types and abundances: 

(1) Mayor et al. (1998), and references therein; (2) Lloyd Evans (1986); (3) Lloyd Evans (1983); (4) Norris & Da Costa (1995); 
(5) Vanture et al. (1994) 



through the strength of the CH G band and the violet CN 
bands. 

The anticorrelation between the CN and CH band 
strengths is equally found in dwarf and giant stars (e.g.. 



Cannon et al. 1998, Briley & Cohen 2001 and references 



therein), and indicates that a C deficiency goes along with 
a N overabundance. Contrary to internal mixing scenar- 
ios, the primordial scenarios generally invoked to account 
for these anomalies extending down to the main sequence 
are difficult to reconcile with C depletions. This difficulty 



may be circumvented if 'C-poor' stars represent in fact 
the primordial C abundance level while C-rich stars have 
accreted carbon, as already suggested in the case of lo Cen 
by Norris & DaCosta (1995). The CN-CH anticorrelation 
should thus be re-interpreted in terms of accretion of ci- 
ther C-poor N-rich matter, or C-rich N-poor matter. In 
the context of the accretion scenario, this situation may 
result depending on whether or not hot-bottom-burning 
(e.g. Boothroyd et al. 1993) was operating in the polluting 
AGB star, as already considered in the simple model de- 
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veloped by Norris & DaCosta (1995). This interpretation 
of the CN-CH anticorrelation requires that ejecta from 
different AGB stars (being either C-poor N-rich, or C-rich 
N-poor) do mot mix in the intracluster gas, a constraint 
difficult to reconcile however with the model of the central 
gas reservoir developed in Sect. ||. 

4.2.2. Na anomalies in main sequence and subgiant 
stars 

Enhancements of Na and Al are common in GC giants, 
and are correlated with N enhancements, and C and O 



depi ctions (|Cottrcll fc Da Costa 1981| , [Norris et al. 1981 



and Da Costa 1997 , and references therein). The first in- 
dication that Na enhancements are observed as well in 
main sequence stars has been provided by Briley et al. 
(1995) for 47 Tuc. The reported enhancements ([Na/Fe] 
~ 0.3 to 0.4 dex) in three CN-strong main sequence 



stars (Afv 1-5) of 47 Tuc (Cottrell & Da Costa 1981 


Gratton et al. 1986; 


Brown & Wallcrstein 1992 


)• 



Very recent observations of turn-off and subgiant stars 
in GCs show that they exhibit clear 0-Na and Mg-Al anti- 
correlations (Gratton et al. 2001), a behavior which had 
previously been seen only in GC giant stars. This result 
cannot be explained by deep mixing scenarios and requires 
some other mechanism, such as accretion of polluted intr- 



acluster gas (Ventura ct al. 2001). 

These anomalies that clearly point against an internal 
mixing scenario do not seem restricted to massive, metal- 
rich clusters like 47 Tuc, since subgiant stars observed by 
King (1998) in the metal-poor cluster M92 ([Fe/H] = - 
2.5) also reveal large Na overabundances (0.8 dex) corre- 
lated with Mg deficiencies. Interestingly enough, these Na 
overabundances are accompanied by modest Ba overabun- 
dances (0.4 dex). 

There are nevertheless indications that qualitatively 
different processes may be at work in metal-poor and 
metal-rich clusters. For example, the constancy of the ra- 
tio of CN-weak/CN-rich stars from the main sequence 
to the RGB observed in 47 Tuc, along with the con- 
stancy of the [Na/Fe] abundances, contrast with the pro- 
gressive Na enhancements and C depletions with increas- 
ing luminosity reported for the more metal-poor clus- 
t crs M13 (iPilachowski et al. 1996[ ), M4 and NGC 6752 
( ^untzcff fc Smith 1991 ). To complicate the picture fur- 
ther, the opposite trend is observed i n M92, since the 
Na overabundances in the subgiants ( King et al. 1998 ) 
are larger than in the RGB-tip giants ( ^hetrone 1996 ). It 
should be remarked that, in the context of the accretion 
scenario, a more relevant distinctive parameter might be 
the cluster mass (47 Tuc is indeed very massive) rather 
than its metallicity. 

In addition, Armosky et al. (1994) has observed that 
in some clusters, the abundances of Y, Ba, Ce, and Nd in 
giant stars are nearly constant as Na/Fe varies widely. 



The different observed behaviors in the Na abundances 
clearly reflects the diversity of the nucleosynthetic pro- 
cesses which produce that chemical element. Sodium and 
s-elements are efficiently produced through proton mix- 
ing by the NeNa cycle (see for example figures 1 and 7 
in Goriely and Mowlavi 2000| ). On the other hand, the 
Na-production during the third dredge-up is not corre- 
lated to the s-process (Mowlavi 1999). Finally, Na is pro- 
duced in massive stars through the combustion of carbon 



(Woosley fc Weaver 1995) 



4.3. Period variations of RR Lyrae in globular clusters 

Cox (1998) has suggested that the slight period decrease 
(2.5 10-5 d over 80 y) observed for the RR Lyr V53 in 
M15, which cannot be attributed to secular evolutionary 
effects, is caused by the accretion of about 10" over k, 
100 y. The corresponding accretion rate (10"^ M© y~^) is 
however much larger than that predicted to result from the 
accretion of gas from a central reservoir in M15. Table |^ 
predicts that a 1 M© star accretes 0.07 M© over the cluster 
lifetime of ~ 10^" y, resulting in an average accretion rate 
of 7 10-12 M0 y-i. 



5. The accretion scenario 

It is worthwhile to attempt to understand the fate 
of all the interstellar matter ejected by massive and 
intermediate-mass stars during the early phases of a GCs 
evolution. On the one hand, it can bring an answer to 
the apparent lack of ISM in present-day GCs and, on the 
other hand, if the efhciency of the accretion process can 
be related to some properties of GCs, it may be expected 
to unravel why GCs exhibit different abundance patterns 
in their main-sequence, turn-off and subgiant stars. 

5.1. Bondi accretion rate 



Following Bondi's ([Bondi 1952[ ) model for spherical accre- 
tion, the rate of mass accretion by a star of mass Ms is 
given by 



dM, 
dt 



-3/2 



(1) 



where A is a constant of order unity, pg is the unperturbed 
gas density, Cs is the sound speed in the gas, and Wrei is 
the relative velocity of the star with respect to the gas. 
In what follows, all masses are given in units of 1 Mq. 

5.2. Mass loss by AGB stars 

First we look at the amount of gas ejected by the AGB 
stars. We assume that the stars eject their mass instan- 
taneously when reaching the AGB phase. The rate of gas 
production by AGB stars is given by 



drrig 



-m, 



dN 



(2) 
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where TOoj is the mass ejected by AGB stars, and dN/dt 
is the rate at which stars reach the AGB phase. We can 
write 



diV 
'df 



dN dm 
dm dt 



(3) 



where dN/dm is the initial mass function (IMF), giv- 
ing the number N of stars born in a given mass inter- 
val {m,m + dm), and dm/dt is the rate at which stars in 
this mass range reach the AGB phase. We choose a simple 
power-law initial mass function: 



dA^ 

- — = Al mci.o m 
dm 



(4) 



where TOci^o is the total initial cluster mass, a is the power- 
law index, and the normalization constant 



ATi = (a - 2)(m 



(5) 



where mi and m^ are the lower and upper mass limits 
of the mass spectrum. For a Salpeter-like IMF, we have 
a = 2.35, whereas if we want to favor the higher mass 
stars we choose a flatter spectrum with, e.g., a = 1.5. The 
evolution time ^ms (in y) of a star of mass m (between 1 
and 10 Mq) along the main sequence can be approximated 



by ( iBahcall fc Piran 1983|) 



log ^MS = 10 - 3.6 log m + (log i 



(6) 



Here we neglect the additional post-main-sequence evo- 
lution time, which is small compared to ^ms- From this 
relation we get 



dm 
dt 



1 



MS 



tMS 2 logm — 3.6 



(7) 



Finally, we adopt the updated results of Weidemann 
(2000) for the initial-to-final mass relation. We fit their 
results: 



0.053 TOi + 0.497 ifl.0<mj<2.5 
0.113 mi + 0.347 if 2.5 < m, < 4.0 
0.073 TO, + 0.507 if 4.0 < m, < 8.0 



(8) 



mi are the final and initial masses of the 
TOcj = mi — mf. Putting together equa- 



where mf and 
star, so that 

tions (^)-(||), we can calculate the rate at which the gas is 
ejected by the AGB stars as a function of time. The total 
amount of gas which has been ejected by the AGB stars in 
the cluster after a given time t is obtained by integrating 
the gas production rate: 



dmg 
dt' 



dt'. 



(9) 



where ti is the time at which the most massive stars (of 
mass TO„) reach the AGB, i.e., logii = 10. — 3.61ogTO,u + 
(logTOu)^. In Figure ^, we show (thin solid line) mg(t) for 
m„ = 8, TO/ = 0.1, and a = 2.35. To illustrate the sen- 
sitivity on the IMF parameters, we also show the results 
obtained for other values of these parameters. Increasing 
the upper limit of the mass spectrum of course increases 



the number of high mass stars, resulting in more matter 
being ejected. Lowering the lower mass limit mi increases 
the total mass situated in low-mass stars, thus decreasing 
the relative number of high-mass stars. By far the most 
sensitive parameter is a, the power-law index of the IMF. 
Choosing a flatter mass spectrum (a = 1.5) considerably 
increases the number of high-mass stars in the cluster, 
leading to a much larger amount of gas ejected into the 
cluster's ISM. With this value of a we find that about 60% 
of the cluster's initial mass can be returned to the cluster 
as gas in lO^^y, whereas for a Salpeter's IMF, this num- 
ber is 20%. We also see that the rate of gas production in 
the cluster is very high at early times, and decreases very 
quickly as the most massive stars leave the AGB phase. 
Note that in this figure mg is the total amount of matter 
which has been lost by the cluster's AGB stars. This is 
different from the amount of gas available for accretion, 
since some of it could have already been accreted, or could 
escape before being accreted. 




Fig. 2. Mass of gas ejected by the AGB stars, normalized 
to the cluster's initial mass. We have used a = 2.35, to; = 
0.1, and TO„ = 8 (thin solid line); a = 1.5, to; = 0.1, 
and TOti = 8 (thin dotted line); a — 2.35, to; = 0.01, and 
to„ — 8 (thick solid line); a = 2.35, to; = 0.1, and to„ = 7 
(thick dotted line). 



If we assume that all the gas ejected by cluster stars 
remains in the cluster and that the only gas-removal mech- 
anism is by tidal sweeping when the GC crosses the galac- 
tic disk, we get a strict upper limit on the amount of gas 
which can be found at any time in the cluster. Even though 
there could have been a very large amount of gas in the 
cluster before its first passage through the galactic disk, 
this amount becomes much smaller afterwards. 
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5.3. Accretion by cluster stars 

The rate of mass accretion by a star of mass Mg is given 
by equation (|l|), where the gas density pg is determined 
by the amount of gas ejected by the AGB stars, but also 
by the ability of the cluster to retain this gas, and by 
the amount of gas which is accreted. We can write this 
equation as 



drris 
"dT 



2.376 10-'' Xmlpv 



10 



(10) 



where m, = M./Mq, r = t/lO^y, vio = (w^^i + 
c2)i/2/(io km/s) and p = /9g/(lMo/pc3). 

The "average" gas density in the cluster can be written 

as 



{rrig 



(11) 



47rr3/3 

where the radius of the gas reservoir is given in pc, 
is the total mass of accreted gas, mesc is the mass of gas es- 
caping from the cluster, either through a continuous wind 
or by tidal forces (for example when the cluster crosses 
the galactic disk). The actual gas density is of course a 
function of r, and we expect it to adopt some smooth pro- 
file peaking at the cluster's center and decreasing rapidly. 
Here we assume that the gas sinks to the cluster center, 
forming a central reservoir of homogeneous density and of 
radius comparable to the core radius re- Typical values for 
Tc range from 0.1 to 10 pc. If we assume that all the clus- 
ter stars have stochastic orbits, the fraction of stars in the 
cluster core is also the fraction of time spent by each star 
in the cluster core. On the other hand, stars that have sunk 
to the center could have smaller velocities and would tend 
to remain in the core for longer periods of time, accreting 
most of the gas (Faulkner & Coleman 1984). Here, we use 
the first hypothesis, keeping in mind that some stars find- 
ing themselves preferably in the core could accrete more 
mass than calculated here, while others would accrete less. 
Using the model of King (1962), the ratio 7 of the num- 
ber of stars in the core to the total number of stars in 
the cluster ranges between 0.3 and 0.1 approximately, for 
concentration parameters c = log(r(/rc) between 1 and 2 
(rt is the tidal radius). We ignore mass segregation effects. 

At a given time, a star is accreting gas at a rate which 
depends on its position in the cluster, its velocity with 
respect to the gas, and its own mass (which itself can vary 
with time as the star accretes or ejects gas). To get the 
total amount of gas which has been accreted by a given 
star of mass at a given time r, we must integrate the 
mass accretion rate over time: 



^a,s(T) = 7 



^ dm . 



dr' 



dr' 



(12) 



where ri = fi/lO^y. We note that in equation (|T^, at any 
time r we have ms(r) — rrisfl + rua^siT) where m^^o is the 
star's initial mass. Finally, the total mass of gas accreted 



by the stars is obtained by integrating m^ s over the mass 
spectrum: 



mair) 



dN ^ 
' dm 



-dm. 



(13) 



Here, for the sake of simplicity, we have made the assump- 
tion that stars do not accrete matter after reaching the 
turn-off point. Indeed, those stars will instead eject mat- 
ter through stellar winds. In equation (13), m^^s and the 
turn-off mass rriTo both depend on time. This is also the 
case for the mass spectrum dN/dm. Indeed, higher-mass 
stars die leaving remnants of lower masses, and some stars, 
especially the lower-mass stars, escape from the cluster by 
evaporation or through external tidal forces. The accretion 
process can also modify the mass spectrum if some stars 
accrete an appreciable fraction of their own initial mass. 
For simplicity, we assume that the mass spectrum remains 
constant with time. Since most of the gas is ejected at early 
times, as shown in Fig (|^), most of the accretion also takes 
place at those early times and we expect that during the 
time interval of highest accretion rate, the variation in the 
mass spectrum is not significant. 

6. Results 



M5 




10 12 14 



Fig. 3. Results obtained for the cluster M5: rua (thin solid 
line), mosc (dotted line) and rrires (thick solid line) as a 
function of time. 

We solve the system of equations (1)-(13). We as- 
sume that the gas forms a homogeneous central reser- 
voir of radius r„ = 



and that its temperature T 



5000 K ( VandenBerg fc Faulkner 1977 ), which gives « 
5.8 km/s. Wc also assume that the stars have chaotic mo- 
tions so that they spend about 20% of their time in the 
cluster core, independently of their mass and average ve- 
locity. The amount of mass accretion is proportional to 



A. Thoul et al.: Accretion of gas by globular cluster stars 



9 



Table 2. Structural parameters for some globular clusters. Unless otherwise noted, the masses M and the velocity 
dispersions CTc are from Mandushev et al. (1991), the central densities pc are from Webbink (1985), the core [tc) and 
tidal (rj) radii are from the catalogue of Harris (1996), the cluster's orbital periods (P) and eccentricities (e) are 
from Dinescu et al. (1999). The escape velocities v^sc are from Madore (1980). and c = log{rt/rc) is the concentration 
parameter (cc stands for "core collapse"). 



NGC 


M 


logPc 




rt 


c 






P 


e 




(IO^Mq) 


(Me/pc») 


(pc) 


(pc) 




(km/s) 


(km/s) 


(io« y) 




104 (47Tuc) 


ir 


3.52" 


0.2" 


48. 


2.40 


10.0 


56.8 


190-193 


0.17 


5139 (w Cen) 


51" 


3.35 


3.7 


64. 


1.24 


16.6 


51.2 


120-123 


0.67 


5272 (M3) 


6.8 


3.86 


1.4 


99.8 


1.85 


4.8 


34.3 


297-321 


0.43 


5904 (M5) 


4.6 


4.11 


0.90" 


56.9" 


1.80 


6.5 


34.9 


722-995 


0.88 


6205 (M13) 


5.8 


3.56 


1.3"= 


47.1" 


1.56 


7.86 


30.4 


429-526 


0.62 


6341 (M92) 


3.9 


4.38 


0.54 


34.9 


1.81 


6.1 


33.1 


201-208 


0.77 


6752 


1.4 


4.47 


0.19 


60. 


2.5(cc) 


4.9 


31.1 


153-156 


0.08 


7078 (M15) 


4.9"= 


6.62'' 


0.64"= 


78.7" 


2.1(cc) 


15.1 


40.9 


242-253 


0.31 



" Murphy et al. 1998 Meylan et al. 1995 " Dull et al. 1997 Gebhardt & Fischer 1995 " Lehman & Scholz 1997. 



Table 3. Results for the clusters listed in Table 2, at t = 14 Gy: ma, total amount of gas that has been accreted by 
the cluster stars, mosc, amount of gas which has been lost from the cluster (swept when crossing the galactic plane), 
TTircs, the amount of residual gas in the cluster (varies between just after crossing the galactic plane and its maximum 



value just before crossing the plane), mknapp the upper limit of the observed dust and ionized gas (Knapp et al. 1995 
1996), and the fraction of mass accreted by a 1 Mq star (last column). 



NGC 


ma/md.o 


mcsc/?Ti,cl,0 


nia/mg 






ma/m,, (mi = 1 Mo) 




(%) 


(%) 


(%) 


(Mo) 


(Mo) 


(%) 


104 (47Tuc) 


21 


0.2 


98 


<7 


1 


80 


5139 (L^Cen) 


0.8 


21 


0.4 


<542 


14 


0.1 


5272 (M3) 


5 


16 


24 


<145 


14 


11 


5904 (M5) 


14 


8 


63 


<140 


7 


40 


6205 (M13) 


4 


17 


18 


<233 


7 


8 


6341 (M92) 


14 


8 


63 


<33 


34 


37 


6752 


20 


1 


95 


<1 


5 


76 


7078 (M15) 


4 


18 


16 


<69 


298 
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the inverse cube of vio and Tc, so that slight variations 
in these parameters induce large differences in the ac- 
creted mass. For the relative stellar velocities, we use the 
"average" value, given by the velocity dispersion in the 
cluster core, CTc- We use the present values for the struc- 
tural and dynamical parameters of the GCs, as listed in 
Table 2, even though those parameters could have changed 
as the clusters evolved chemically and dynamically. These 
clusters have been selected on the basis of their high 
escape velocities (these clusters are tightly bound and 
therefore more likely to retain gas) and/or because abun- 
dance anomalies have been observed for stars belonging to 
these clusters. The reservoir of gas is emptied each time 
the GC crosses the galactic disk, i.e. every P/2 where 
P is the cluster's orbital period. This is the only con- 
tribution considered here for the term mosc appearing in 
Eq. |ll|. Other removal processes, like type I supernovae, 
kinetic effects of strong winds from millisecond pulsars 
( ^pergel 199"l| , [Freire et al. 2001 or a photoionized clus- 
ter wind driven by the UV radiation from hot stars (see 
Sect. I and [Knapp et al. 1996| , [Freire et al. 2001| ), were 
not included in mcsc. Any such process would decrease 



the efficiency of the present accretion scenario. The re- 
sults presented here must therefore be considered as upper 
bounds. 

Finally, we use a — 2.35, m; — 0.1, m„ = 8, and A = 1 
in all cases. The results at i = 14 Gy are listed in Table 3, 
where we give the total amount of gas which has been 
accreted by the stars, m^, the amount of gas lost from the 
cluster (swept when crossing the galactic plane), rriesc, the 
maximum amount of residual gas in the cluster, rrires, and 
the fraction of mass accreted by a 1 star. 

Over the lifetime of the cluster, the total amount 
of matter injected into the GC-interstellar medium by 
the AGB stars is about 20% of the cluster's initial mass 
(see Fig. 2). In the most concentrated clusters, most of 
this gas is accreted by the cluster's lower-mass stars. 
Indeed, in 47Tuc, NGC6752, M5, and M92, more than 
60% of the gas ejected by the AGB stars is accreted 
by the cluster stars. In those clusters, 1 M0 stars can 
accrete an appreciable fraction of their initial mass. The 
envelopes of those stars will reflect the composition of 
the intracluster medium rather than the composition of 
their interior, even if a physical mechanism is at work 
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to induce mixing with the deeper layers of the star. We 
note here that these clusters are rich in observed stellar 
composition anomalies (see Sect. 4). The only cluster for 
which we get a very low amount of accretion is uj Cen. 
This is due to a combination of factors: w Cen has a low 
concentration parameter compared to the other clusters 
we have selected. It also has a rather large core radius, 
and a high core velocity dispersion, as well as frequent 
disk crossings. On the other hand, lo Cen is a peculiar 
GC in many other respects and should probably not 
be considered as a test bed for the accretion scenario, 
as many studies now view to Cen as a conglomerate of 
different subsys tems or even a s the nucleus of a dwarf 



galaxy (e.g.. 


Lee et al. 1999|, 


Hilker & Richtler 2000 


Majewski et al. 200C, 


Pancino et al. 2000 


)■ 



The upper limits we get for the current residual intra- 
cluster gas are usually larger than the observational limits 
determined by Knapp et al. (1995, 1996), but of course 
the exact amount of gas predicted in the cluster will de- 
pend on the time elapsed since its last passage through the 
galactic plane. Moreover, other gas sweeping mechanisms 
that were not included in the present model are possibly 
at work, as discussed above. 

The results for m^, rriesc and rrires as a function of 
time are shown in Fig. ^ for one particular cluster, M5. As 
expected, we find that most of the accretion takes place at 
early times, when the amount of gas ejected by the AGB 
stars is the largest. The amount of residual intracluster 
gas increases with time in between each passage through 
the galactic plane, but becomes very small at late times. 

7. Conclusions 

It was shown in this paper that accretion by low-mass 
GC stars of the gas ejected in the intracluster medium by 
moderately massive stars may be quite efficient. If enough 
mass is accreted, it can lead to major alteration of the 
stellar surface composition. This supports the EASE sce- 
nario and provides a plausible explanation for the lack of 
intracluster gas, and for some of the abundance anomalies 
observed in GC stars. 
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